AD-A164 


TECHNICAL  REPORT 


A  HIGH-ORDER  LAMINATED  PLATE  THEORY 
WITH  IMPROVED  IN-PLANE  RESPONSES 


by 


A.  Toledano^^  and  H.  Murakami^^ 

Department  of  Applied  Mechanics  and  Engineering  Sciences 
University  of  California  at  San  Diego 
La  Jolla,  California  92093 

DECEMBER  1985 

Prepared  for 

Office  of  Naval  Research 
800  North  Quincy  Street 

Arlington,  Virginia  22217 


This  document  has  been  approved 
fot  public  release  and  sale;  its 
distribution  is  unlimited. 


Research  was  supported  by  the  Office  of  Naval  Research  under  (rant  N000I4.84-K-O468  to  the  University  of  California  at  San 
ego. 

Research  Assistant. 

Assistant  Professor  of  Applied  Mechanics. 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


«.  TlTwC  (an«  SukiKU; 


A  High-Order  Laminated  Plate  Theory  with 
Improved  In-Plane  Responses  1) 


CATkwOC  NUMttCN 


PERroRMiMO  e«».  nK»OAT  NuMacA 


7.  AjTrtCA^»> 


A.  Toledano  and  H.  Murakami 


a.  coNiNACt  ON  cramt  NUMacN^ki 


N00014-84-0468 


a.  NtaToNiK'kC  oncakizatioa  kami  and  aocnczs  J5tAV*oVic^iII[iT\'u««eEN/' 

Dept  of  Applied  Mechanics  4  Engineering  Sciences 
University  of  California,  San  Diego 
La  Jolla,  California  92093 


i>.  cosrnoLL'no  ofice  name  and  adoncss 

Mechanics  Division,  Solid  Mechanics  Program 
Office  of  Naval  Research 

800  N.  Quincy  St.,  Arlington,  VA  22217 _ 

*■  MONiTONinC  aCEkCy  name  a  AOCMCiS.'!/ IrwK  CtttttHlna  Ottln)  j  '»  SEClNiTV  CuAii.  f»/ »(* 


IE.  NENONT  DATE 

December  1985 


•  I.  NUMNEN  or  naCES 

41 


rr»a  DECk AtSiTiC ATiON  DONnSNADinC 
•CmEOulE 


ta.  OiSTMlBUTlON  statement  .‘A/I  N*psr.v 

Approved  for  public  release;  distribution  unlimited. 


17.  eiSTNiauTtON  statement  (•!  *Aa  aSAIiAct  tnft»4  tm  OIacA  JO,  M  tlllttml  Nw> 


IS.  NET  «ONSS  (C»ntinut  O.-'  iftrn  a/N*  II  i>acaa»*rr  an«  INantIf/.Sr  WaeA  niaiksr; 


Laminated  Composite  plate,  shear  deformable  plate  theory 


1C.  ASSTNACT  fCani/nua  an  la^araa  alNa  II  nacaaaaTr  miN  INamfly  Or  OlaaA 

In  order  to  improve  the  accuracy  of  the  in-plahe  responses  of  the  shear  deform¬ 
able  laminated  composite  plate  theories,  a  new  high-o^dgr.  theory  was  developed  . 
based  upon  Relssner's  new  mixed ■ Variatiortal  Principal.  To  this  end,  a  zig-zag  , 
shaped  C  Junction  and  Legendre  polynomials  were  introduced  into  the  approximate  t 
in-plane  displacement  distributions  across  the  plate  thickness.  The  accuracy  of  i 
the  present  theory  was  examined  by  applying  it  to  the  cylindrical  bending  problem? 
of  laminated  plates  which  had  been  solved  exactly  by  Pagano(1970) .  A  comparison  | 
with  the  exact  solutions  obtained  for  several  symmetric  and  asymmetric  cross-ply  i 

Irrates  that  LllH  BPSSEBt  tnffflPV  SC curateiy, estimates^iu-plane  responses,  even  for 
1/7T  •n<T>aii  Of  •  MOV  •»  IS  ONEOLETt  small  span-to-thicriness  ratios. 


1.  INTRODUCTION 


The  increasing  use  of  composite  materials  as  thick  laminates,  in  aerospace  engineering  and  in 
automotive  engineering,  has  clearly  demonstrated  the  need  for  the  development  of  new  theories  to 
efficiently  and  accurately  predict  the  behavior  of  such  structural  components.  The  intrinsic  hetero¬ 
geneity  and  anisotropy  of  these  composite  structures  as  evidenced  in  the  stacking  of  several  fibrous 
layers  and  in  the  high  discontinuity  in  material  properties  across  the  interfaces,  make  the  classical 
theories  of  plates  and  shells  inadequate. 

The  inspiration  and  guidelines  for  the  subsequent  attempts  have  stemmed  from  Pagano's  works 
[1,2,31  where  the  exact  elasticity  solutions  for  the  problems  of  cylindrical  bending  and  simply  supported 
rectangular  plates  were  given.  Pagano  showed  the  importance  of  incorporating  the  effect  of  transverse 
shear  deformations  in  order  to  accurately  estimate  the  plate  lateral  deflection  and  the  need  to  improve 
upon  the  thickness  variation  of  the  in-plane  displacements,  which  are  assumed  to  be  C'  linear  functions 
in  both  classical  plate  theory  (CPT)  and  Reissner-Mindlin  plate  theory  (FSD). 

The  first  attempt  to  develop  a  general  linear  laminated  plate  theory  is  credited  to  Yang,  Norris  and 
Stavsky  [4] .  Their  theory  is  an  extension  of  the  Reissner-Mindlin  homogeneous  plate  theory  as  applied 
to  an  arbitrary  number  of  bonded  anisotropic  layers.  Whitney  and  Pagano  [5]  extended  Yang,  Norris 
and  Stavsky's  work.  An  important  conclusion  drawn  from  their  analysis,  which  was  also  emphasized 
later  by  Whitney  (6),  is  that  the  inaccuracies  of  the  classical  plate  theory  at  low  span-to-thickness  ratios 
for  determining  in-plane  stresses  are  not  alleviated  by  the  introduction  of  shear  deformations.  Whitney 
[6]  obtained  in-plane  displacements  by  integrating  the  transverse  shear  strains  deduced  in  [S].  This 
resulted  in  a  higher  order  approximation  which  accurately  predicted  in-plane  strains,  but  the  resulting 
modified  stresses  did  not  necessarily  satisfy  the  original  plate  equilibrium  equations. 

Since  then,  other  high-order  laminated  plate  theories  have  been  proposed  that  account  for 
transverse  shear  strains.  Of  these,  the  Lo,  Christensen  and  Wu  [7]  and  the  Reddy  [8]  high-order 
models  have  served  as  the  foundation  for  the  present  theory.  In  their  paper  [7],  Lo,  Christensen  and 
Wu  used  appropriate  higher  order  terms  in  the  power  series  expansions  of  the  assumed  displacement 
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field  which  was  proposed  by  Hildebrand,  Reissner  and  Thomas  (8]  On  the  other  hand,  Reddy  [9] 
imposed  the  condition  of  vanishing  transverse  shear  strains  on  the  top  and  bottom  surfaces  of  the  plate. 
However,  this  theory  does  not  satisfy  the  continuity  condition  of  ti-ansverse  shear  sticaKS  at  the  inter¬ 
faces. 

The  objective  of  the  present  paper  is  to  improve  the  approximation  of  in-piane  variables  in  lam¬ 
inated  plate  theories.  In-piane  displacements  and  bending  and  stretching  stresses  are  considered  pri¬ 
mary  quantities  in  any  approximate  laminated  plate  analysis:  trs.nsverse  stresses  are  only  of  secondary 
importance  since  they  are  an  order  of  magnitude  smal  er  than  the  primary  bending  and  stretching 
stresses.  By  using  a  new  mixed  variational  principle  proposed  by  Reissner  (10),  the  present  theory  is  a 
high-order  model  which  improves  upon  existing  theories  by  including  in  the  assumed  in-plane  displace¬ 
ment  variations  across  the  plate  thickness;  1)  a  zig-zag  shaped  function  as  detailed  by  Murakami 
(111;  and,  2)  Legendre  polynomials.  The  advantage  of  using  Reissner's  new  mixed  variational  principle 
is  that  it  automatically  yields  the  appropriate  shear  correction  factors  for  the  transverse  shear  constitu¬ 
tive  equations.  Another  attractive  feature  of  the  proposed  theory  is  that  the  number  of  equations  to  be 
solved  is  not  increased  as  the  number  of  layers  becomes  larger  and  larger.  A  comparison  of  the  pro¬ 
posed  theory  with  Pagano's  exact  elasticity  solution  for  symmetric  and  asymmetric  laminated  plates  in 
cylindrical  bending,  shows  that  in-plane  displacements  and  stresses  are  accurately  predicted  by  the  inclu¬ 
sion  of  the  zig-zag  shaped  function  and  the  Legendre  Polynomials. 
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2.  FORMULATION 

Consider  an  AT-layer  laminated  composite  plate,  shown  in  Fig.  I,  with  principal  axes  coinciding 
with  a  Cartesian  coordinate  system  such  that  the  X}-axis  is  perpendicular  to  the  plane 

defined  by  Xi  and  X2.  The  following  notation:  (  1,2,.., N  will  designate  quantities  associated 

with  the  -layer.  The  thickness  of  each  layer  is  where  h  is  the  total  thickness  of  the  plate. 

The  volume  fractions  satisfy  the  relation 

-  1  (1) 

Unless  otherwise  specified,  the  usual  cartesian  indicia)  notation  is  employed  where  latin  and  greek 
indices  range  from  1  to  3  and  I  to  2,  respectively.  Repeated  indices  imply  the  summation  convention 
and  (  is  used  to  denote  partial  differentiation  with  respect  to  x, . 

With  the  help  of  the  foregoing  notation,  the  governing  equations  for  the  displacement  vector  u/*'’ 
and  stress  tensor  cr,}**  associated  with  the  /c**’-layer  are; 

a)  Equilibrium  Equations 

+  (2) 

where  ft  are  the  body  forces; 

b)  Constitutive  Equations  For  Orthotropic  Layers 
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where  C,j  are  the  elastic  constants  and  C.j  Uj“  1,2,6)  represent  the  reduced  stiffnesses  introduced  by 
Whitney  and  Pagano  [51; 


■  •  vj 


'• 
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c)  Strain-Displacement  Relations 


d)  Interface  Continuity  Conditions 

e)  Upper  and  Lower  Surface  Stress  Conditions 


<ri*>-7;+  onxj--^ 


tr^'>  ^  Tr  on  X3  -  -  . 


The  objective  in  developing  a  new  laminated  plate  theory  is  twofold;  first,  to  improve  the  assumed 
variation  of  in-plane  displacements  through  the  thickness  of  the  plate  and  second,  to  include  the  effect 
of  transverse  shear  deformation.  In  order  to  carry  out  this  task,  Reissner's  new  mixed  variational  prin¬ 
ciple  (10]  was  applied  to  the  -layer  composite  plate  whose  middle  surface  occupies  a  domain  D  in  the 
Xj,  Xj-plane; 

j  I  J  (  -OISt jj’ j  dxy  dx,  dx^ 


dxi<ic2+  I 

I  r  Ti  dxy 

io 

r  1  * 

//h 


*'’(x,,X2,-|)r/  -  8w/'^>(xi,X2--j)  rrj  dx,  dx2 


V  V  V 


•> 

v"*!- • 

I." 


where  dDr  denotes  the  boundary  of  D  with  outward  normal  v„  on  which  tractions  T  are  specified  and 
represents  the  xj-domain  occupied  by  the  k*^-layer.  Also  pj,  (...)  implies  the  appropriate  right- 
hand  side  of  (3b).  Due  to  the  nature  of  Reissner's  mi.xed  variational  principle,  Eqs.  (3a)  are  taken  to 
be  the  definitions  of  used  in  connection  with  (7). 
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3.  TRIAL  DISPLACEMENT  FIELD,  TRANSVERSE  AND  NORMAL  STRESSES 

The  high-order  laminated  plate  theory  which  takes  into  account  the  effect  of  transverse  shear 
strains,  is  obtained  by  including  the  Legendre  polynomials  of  order  n—  1,2,3  with  respect  to  the  X3- 
coordinate  to  a  zig-zag  in-plane  displacement  variation  of  amplitude  5, 0r|,X2)  across  the  plate  thickness. 

The  appropriate  trial  functions  used  in  connection  with  Reissner's  mixed  variational  principle 
Eq.  (7)  are  taken  to  be; 

a)  Trial  Displacement  Field 

M,<*>(X|,X2,Xj)  -  +  (■|)'F,Cx,,X2)Fi(()  +  S,  (*i,X2)(-  D* 

(8) 

+  ('j)'f,Cx„X2)/»2({)+  (■j)V,Cx,A2)^3Q:) 

2xi 

where  {  =  and  ({)  are  the  Legendre  Polynomials  of  order  n.  It  is  also  understood  that  =  0. 

is  a  local  X3-coordinate  system  with  its  origin  at  the  center  xj§*  of  the  /t**’-layer,  i.e. 

x/»>  SX3-xi^>  .  (9) 

Eq.  (8)  may  be  regarded  as  a  superposition  of  a  zig-zag  function  and  the  cubic  variation  as  proposed  by 
Lo,  Christensen  and  Wu  [7],  with  the  exception  that  here  Legendre  polynomials  are  used  instead  of 
single  powers  in  X3; 

b)  Trial  Transverse  and  Normal  Stresses 
Ti*’(Xl,X2.*j)  “  (?a*’  0Ci,X2)Fi(z)  +  (xl,X2)f2fe)  + 

(10a) 

+  iri*-”  (x,,x2)  +  rj*>Cx,,x2)iF4(r)+  iri*-'>(x„x2)-  rj*'(x„x2)lF5(‘*)  ; 

TjJ’(Xl,X2A3)  -  Oj*’  Cxi,X2)Fi(z)  +  /?/*'’  (Xi,X2)f«(z)  +  CXirX2)F3(z)  +  /]*' 

(10b) 

+  iT-j*-”  Cx,,X2)  +  ri*’  Cx,,X2)]F4(z)  +  iri*-”  (x„X2)  -  ri*>  Cxi,X2)lfg(r) 

where 


Fs(z)  ~  IQz^  - ^  2  ,  F,(z)  ~  j;^  062^  -  \*z^ ^  z) 


/•7a)--^^^(112r'-4(fe'+3z)  ,  /•,(r)-12fc*-35z^+-^2 

xj**  1  1 

Also  ,  (G<*>  ,/?,<*',  /<*> )  s  f  (1,  xj*>  .  xi*^  )  Ti^>  dxy  (12a] 

/}*>  =  f  Xj<*>'  T,<f>  ic,  (12b; 

><<*> 

In  (10)  and  7/*’  are  the  values  of  T3,  at  the  top  and  bottom  surfaces  of  the  fc***  layer  respec 


tively.  From  (6) 


f/o)  -  7;+  and  -  Tr  • 


The  degree  of  the  polynomials  (z),  1-I-8,  appearing  in  (10a,b)  is  consistent  with  the  order  of  trun¬ 
cation  in  the  assumed  expansions  (8)  for  the  displacement  u/**. 
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4.  LAMINATED  PLATE  EQUATIONS 

Substituting  (8)  and  (10)  into  (7),  using  Gauss’  Tlieorcm  and  the  orthogonality  property  of  the 
Legendre  polynomials  one  obtains; 
a)  Equilibrium  Equations: 

+  04a) 

-  Nj,  + 1  (t;*  +  r-)  +  -  0  (i4b) 

Zai^  -  -  (Ti^  -(-1)'^  7}-)  (14c) 

-  3A/j,  +  (r  -  7]-) +  Ff-~0  (14d) 

FfiaJi  -  (51.3.  +  -T  ^i.)  +  -T  (^-  +  7-a  )  +  /•;-  0  (14e) 


I-  • 


where 


^a0-*  ^a0  ^  ^  ^ 
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I  fi  ,  fi  ,  r,  ,  r,  ,  /-I  j  *«i 


2xj*>  ,h 


tl,  j  Pi(0,(-  D*  ,  (7)'  ,  (-f  [^U)]  05a,b) 


,  Mji  ,Kji  .,  Z}i,Li,)  —  1,  2  7*\(0..(  D*  -  (  0*  ^(*)^  ■’^3***  > 


(■j)'/’2((:)  rif'dxj  ; 


b)  Constitutive  Equations: 

•  For  Transverse  Stresses 


(T-j*-"  +  rj*’ )  -  I  A  ej*’  [  i/3,  +  'P,  +  5.(-  D* 
+  h  +  3f  J  +  -^  (3  -  -J-  +  -^  (5  «  |  )<A„ 


.V 

A  A  A 
A’A*>' 
s.'  v' 

A» 


(16b) 
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n'*»  Cj*’ (^3.  +  5^J 


(16c) 
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(16d) 
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>(*) 


'(*+1) 


•rik) _ 2. 

•  n  ^ 


,(*+1) 


Ut-U 


r(*+i» 

*  tt 


•  For  Normal  Stresses 


ij*)  - 


(/|<‘'/|)Z  30 
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where  in  (16a,b,c)  and  (17a,b,c,d)  k  ranges  from  1  to  N  while  in  (16d)  and  (17e)  k  ranges  from  1  to 
W-1) .  Also,  no  summation  on  a  is  implied  in  (16)  and 


=  S„i  +  8„2  Cif»  ;  n  =  x^^/h 


(18) 


u 

U7.2 

'<'2.2 

s 

- 

■5i.i 

■52.2 

1 

f2.2 

<Ai.i 

4>2.2 

(19) 


By  solving  (16)  and  (17),  (?/**,  7?/*’,  7/*’,  /j*’  and  f/*’  are  obtained  in  terms  of  C/,, and 
and  their  derivatives.  As  a  result,  the  quantities  Na,  Mi,,  Ka,  La  of  Eq.  (15c)  can  be  determined 
as  functions  of  these  displacement  variables.  Such  expressions  will  automatically  include  the  appropri¬ 
ate  shear  correction  factors  by  virtue  of  the  Reissner  mixed  variational  principle. 


The  equilibrium  equations  (14)  are  supplemented  with  the  following  suitable  boundary  conditions: 


specify 

u. 

or 

^at  » 

(20a) 

specify 

or 

^at  ' 

(20b) 

specify 

s, 

or 

^ai  ' 

(20c) 

specify 

or 

^at  » 

(20d) 

specify 

<t>a 

or 

(20e) 

The  remaining  constitutive  equations  for  Map,  Zap,  Lap  and  Pap  are  obtained  by  substituting 


(3a),  (4),  (8)  and  (10b)  into  (15a)  to  yield: 


if 
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symmetric  [I^I 
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(21) 


where  DVn  )V22  N12V,  It/i.i  Ui^Uxi^  analogous  expressions  for  'I',  ...,  P, 

^  [P^l  are  3  X  3  matrices,  (Cl^*^  is  a  15  x  5  matrix  an  ...,  V*  are  1  x  4  vectors,  which 


are  given  in  the  Appendix. 
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5.  CYLINDRICAL  BENDING  OF  LAMINATED  PLATES 

In  order  to  test  the  accuracy  of  the  present  theory,  cylindrical  bending  of  composite  plates  under 
sinusoidal  loading  is  considered.  The  plate  is  simply  supported  at  the  ends  jti  >  0  and  /  and  is  infinitely 
long  in  the  xj-direction.  The  prescribed  boundary  conditions  on  the  top  and  bottom  surfaces  of  the 
plate  are; 

rf  *  0 , 7^  -  9  sin  on  xj  -  -j  (22a) 

7T-7T-0  onx,---|.  (22b) 

The  boundary  conditions  for  the  simply  supported  ends  are,  from  (20); 


at  xj  —  0  ,  / 

(23a) 

iVii  -  M\\  -  Z\\ -  Ill  -  /’ll  “  0 

at  xi  «  0  ,  /  . 

(23b) 

Using  surface  boundary  conditions  (22),  the  equilibrium  equations  (14)  for  cylindrical  bending  reduce 
to; 


o 

1 

(24a) 

'^13,1  +  9  sm  — ^  -  0 

(24b) 

Mi.i  “  'Y3,  -  0 

(24c) 

M\y\  -  iV33  +  <7  sin  — ^  —  0 

(24d) 

0 

1 

1 

N 

(24e) 

Zi3,i  -  ^33  ~  9  sin  -  0 

(240 

1 

1 

0 

(24g) 

From  the  boundary  condition  -•  0  at  X|  —  0,/,  Eq.  (24a)  implies  that 


/Vii-O  . 

(25) 

Next  Eqs.  (15a,c)  are  expressed  in  terms  of  the  displacement  variables  t/j . fj.  To  this  end,  the 

constitutive  equations  (16)  and  (17),  for  the  cylindrical  bending  analysis,  can  be  rewritten  in  the  fol¬ 
lowing  vector  form: 

Qi~ 

(26a) 

AlB,Jr,-X2 

(26b) 

(26c) 

in?.  1  Gi  +  {  17^1 1^1  - 1 ITlPj I7i  -  A ic, )  r, 

(26d) 

and 

~  73  ^  I^  1 1  r3 "  ?i 

(27a) 

7(l3--^73-7^lfl.lr3-?2 

(27b) 

Qi~  ■j^73+'|’*l^llr3“?3 

(27c) 

I  ?■- if 

(27d) 

in?3lG3+ jl7K3M3-|-;jT  ln?3l73--f|-;Jr  itr,)/,- 

/ilCjlfj  (27e) 

8//  -  32 

/-1,3  and/,=  -^  /,  . 


where  7,  = 


(28) 


The  matrices  lAi) .  IC3I  and  vectors  X, . 54  are  given  in  the  Appendix.  The  vector  equations 

(26a,b,c)  and  (27a,b,c,d)  have  ^-components,  while  the  vector  equations  (26d)  artd  (27e)  have 
W-  1)  components.  Matrices  [A\],  ■■■,  (C3I  depend  on  the  volume  fractions  end  elastic  con¬ 
stants  Cif ,  C/jf’  and  Cjf  * ,  while  the  vectors  ^4  contain  the  displacement  variables  Ui ^3 . 

Eqs.  (26)  are  easily  solved  by  substituting  Qi,  and—  7,  in  terms  of  T,  from  (26a,b,c)  into 
(26d).  This  yields  a  new  equation  involving  Tj  only,  which  can  thus  be  solved  for  Ti.  Then  by  back 
substitution  expressions  for  Q],  R\  and  J\  in  terms  of  Xj,  Xj  and  X3  are  obtained.  Proceeding  in 
a  similar  manner  with  (27a,b,c,d)  Qi,  /?3,  Jy  and  ly  in  terms  of  *t|,  kj,  K3  and  K^  are  deter¬ 
mined.  These  expressions  are: 


[(jl/l-HO,!)  -(-jI/1-2L40i1)  [ki| 


4  «1-  IfiOl  1(^1-  2X3)-!-  (l/l-  l^/iiDXj 


(y  {/)  -  UQy])  -  (yl/J  -  2UQy])]  k 


{|l/)-2U03l)  -  (yl/l-4U(?3l) 


,  lARy] 
ley'll  2\ARy\ 


±J 

* 

^^U]+-^lBRy])\  PI 


1  101BO311 

IT  16U?0j1 


where  [/)  is  the  x  identity  matrix  and 


[AQ,]  IARi]  l4 il  I  ] 

[BOA  “  (flil  [77?,]|  /-1.3 
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with 


(7r,l-  (41710, 1L4,1  +  (r/i,llB,l  +  IC,])-'  (31b) 

I7y,l  -  (4(7103lU,l  -  (77? 3)15,1  +  (C,l)-'  (31c) 


\AQi  1 .  IBR,  ]  are  ,V  x  ,V  matrices,  while  (TT) )  are  (,V—  1)  x  CV—  1)  matrices.  By  inserting  (29) 

and  (30)  into  (ISc)  and  (21)  the  appropriate  constitutive  relations  for  the  cylindrical  bending  problem 
in  terms  of  the  displacement  variables  7/,, and  their  derivatives  with  respect  to  jt,  are  obtained. 

The  form  of  the  dependence  on  the  displacement  variables  f/,,  ...,  (3  of  the  constitutive  equations 
thus  obtained  and  the  nature  of  the  applied  load  suggest  the  following  expressions  for  the  displace¬ 
ments: 


Ui 

/li/l 

Ui 

hUi 

% 

X, 

^^3 

Si 

•• 

AS, 

cos  rr  —  and 

S3 

* 

ASi 

(l 

Uh 

6 

yft 

1 

(32) 


where  the  quantities  are  nondimensional  by  definition.  It  is  easily  proven  that  the  boundary  condi¬ 
tions  (23)  are  satisfied  when  (32)  are  substituted  therein. 

Finally,  inserting  (32)  into  the  constitutive  equations  obtained  in  the  manner  described  above  and 
these  in  turn  into  the  equilibrium  equations  (24)  and  (25)  yields  a  system  of  nine  algebraic  equations 
with  the  nine  nondimensional  quantities  t/,,  I3  as  unknowns.  This  system  is  conveniently  written 
in  matrix  form  as 


[B]U^F 

(33) 

where  f/- (£/, 'P,  S,  <i,  f/3  ^3  5313)^ 

(34a) 

F-  \0,g,0.^q,0,-q,0,^q,0V 

(34b) 

and  [5]  is  a  9  X  9  matrix. 


6.  NUMERICAL  RESULTS 


-  18- 


represented  by  a  dashed-dotted  line  and  dotted  solid  line,  respectively.  Symmetric  3,  S  and  9-ply  lam¬ 
inates  and  asymmetric  4  and  8-ply  laminates  were  examined,  to  test  the  present  theory. 

For  a  symmetric  3-ply  laminate  (0/90/0)  with  layers  of  equal  thickness.  Table  1  shows  the  values 
of  the  central  deflection  U3  obtained  from  the  different  theories  for  a  span-to-thickness  ratio  5  of  4  and 
6.  As  observed  the  present  high-order  theory  correctly  predicts  the  central  deflection  U3  to  the  first  two 
decimal  digits,  while  the  first  order  zig-zag  model  gives  a  better  result  than  LCW.  The  variation  of  the 
in-plane  displacement  Uj  across  the  plate  thickness  is  compared  in  Fig.  2a  for  5  —  4,  where  it  is  seen 
that  the  curves  for  the  present  theory  and  the  exact  solution  are  almost  identical.  This  improvement  is 
also  reflected  in  the  variation  of  the  in-plane  stress  an  across  the  plate  thickness,  as  shown  in  Fig.  2b. 
Very  close  agreement  is  found  between  Pagano's  exact  solution  and  the  present  theory,  which  has 
improved  upon  Lo,  Christensen  and  Wu's  high-order  theory,  especially  at  and  in  the  neighborhood  of 
the  interfaces. 

The  present  theory  was  next  tested  for  a  symmetric  5-ply  laminate  (0/90/0/90/0)  with  layers  of 
equal  thickness.  The  central  deflection  U3  for  span-to-thickness  ratio  S  of  4  and  6,  is  shown  in 
Table  1  where  close  agreement  with  the  exact  solution  is  observed.  The  variations  across  the  plate 
thickness  of  in-plane  variables  u{**  and  a/f’  are  compared  in  Figs.  3  and  4.  The  curves  for  the  present 
high-order  theory  and  the  exact  solution  are  again  almost  identical.  In  particular,  it  is  seen  that  the 
present  theory  has  considerably  improved  upon  Lo,  Christensen  and  Wu's  model  in  the  interior  layers 
of  the  plate. 

To  further  assess  the  accuracy  of  the  present  high-order  theory  the  more  diflicuit  case  of  a  sym¬ 
metric  9-layer  cross-ply  laminate  (0/90/0/90/0/90/0/90/0)  was  considered.  The  0®  layers  have  equal 
thickness  h/10  while  the  90®  layers  have  equal  thickness  h/S.  The  results  for  the  central  deflection  C3 
are  given  in  Table  1  for  5  4  and  6  where  again  close  agreement  with  the  exact  solution  is  observed. 

The  variations  across  the  plate  thickness  of  the  in-plane  displacement  U|  and  normal  stress  are 
shown  in  Figs.  5  and  6,  for  S  ■■4  and  6  respectively.  There  the  discrepancies  between  the  first  order 
zig-zag  theory  and  the  exact  solution  are  more  pronounced  than  in  the  3-  and  S-layer  cases,  as 


r-- 


expected.  However,  the  results  of  the  present  theory  are  still  very  good  when  compared  to  the  exact 
solution. 

Finally,  asymmetric  4  and  8  cross-ply  laminates,  with  layers  of  equal  thicknen,  were  examined. 
The  present  theory  predicts  accurately  the  central  deflection  uj.  These  results  are  given  in  Table  2  for 
span-to-thickness  ratio  5  of  4  and  6.  The  variation  across  the  plate  thickness  of  the  in-plane  displace¬ 
ment  u  f**  and  normal  stress  5/*’  are  shown  in  Figs.  7,  8  and  9  for  S  *  4  and  6.  From  the  curves  for 
u/''\  it  is  seen  that  the  first-order  zig-zag  theory  deviates  significantly  from  the  exact  solution  at  the 
bottom  layer  of  the  plate.  On  the  other  hand,  the  discrepancies  between  LOW  and  the  exact  solution, 
for  both  u  I*  *  and  o  If  ’  are  more  pronounced  in  the  interior  layers  of  the  plate,  while  the  present  high- 
order  theory  is  still  in  very  good  agreement  with  the  exact  solution. 


7.  CONCLUSION 


A  high-order  laminated  plate  theory,  which  accurately  predicts  in-plane  responses  of  symmetric 
and  asymmetric  laminates,  was  developed  with  the  help  of  Reissner's  new  mixed  variational  principle 
[10].  The  improvement  was  achieved  by  including  a  zig-zag  shaped  C’‘  function  in  the  in-plane  dis¬ 
placement  variations  across  the  plate  thickness,  as  proposed  by  Murakami  III],  while  the  non-linear 
variation  is  accounted  for  by  using  Legendre  Polynomials.  The  accuracy  of  the  theory  was  examined 
for  the  case  of  cylindrical  bending  of  an  infinitely  long  strip  and  compared  with  the  exact  elasticity  solu¬ 
tion  given  by  Pagano  [1].  The  results  for  the  central  deflection  and  in-plane  displacements  and  normal 
stresses  for  several  symmetric  and  asymmetric  cross-ply  laminates  indicate  that  the  theory  very  accu¬ 
rately  predicts  these  in-plane  responses  even  for  small  span-to-thickness  ratios.  In  all  the  cases  con¬ 
sidered,  the  proposed  theory  gave  better  in-plane  responses  than  the  Lo,  Christensen  and  Wu  high- 
order  theory,  especially  in  the  interior  layers  of  the  plate.  It  was  also  observed  that  for  symmetric  lam¬ 
inates,  the  first  order  zig-zag  model  [II]  predicts  more  accurately  the  central  deflection  than  the  Lo, 
Christensen  and  Wu  high-order  theory. 


ACKNOWLEDGEMENT 

The  research  was  supported  by  the  Office  of  Naval  Research  under  Grant  N00014-84-K-0468  to 
the  University  of  California  at  San  Diego.  The  authors  are  grateful  to  Dr.  A.  Kushner  for  his  support. 
Also,  the  author  is  thankful  to  Professors  E.  Reissner  and  G.A.  Hegemier  for  their  advice  and  com¬ 
ments  on  the  manuscript. 


REFERENCES 


1.  Pagano,  N.  J.,  ’Exact  Solutions  for  Composite  Laminates  in  Cylindrical  Bending,*  J.  Comp.  Mat., 
i.  398-411  (1970). 

2.  Pagano,  N.  J.,  ’Exact  Solutions  for  Rectangular  Bidirectional  Composites  and  Sandwich  Plates,’  J. 
Comp.  Mat.,  4,  20-34  (1970). 

3.  Pagano,  N.  J.  and  Hatfield,  S.J.,  ’Elastic  Behavior  of  Multilayered  Bidirectional  Composites," 
AIAA  Journal,  10,  931-933  (1972). 

4.  Yang,  P.  C.,  Norris,  C.  H.,  and  Stavsky,  Y.,  ’Elastic  Wave  Propagation  in  Heterogeneous  Plates,’ 
Int.  J.  Solids  Struct.,  2,  665-684  (1966). 

5.  Whitney,  J.  M.  and  Pagano,  N.  J.,  ’Shear  Deformation  in  Heterogeneous  Anisotropic  Plates,” 
ASME  J.  Appl.  Mech.,  37,  1031-1036  (1970). 

6.  Whitney,  J.  M.,  "Stress  Analysis  of  Thick  Laminated  Composite  and  Sandwich  Plates,’  J.  Comp. 
Mat.,  6,  426-440  (1972). 

7.  Lo,  K.  H.,  Christensen,  R.  M.,  and  Wu,  E.  M.,  ’A  High-Order  Theory  of  Plate  Deformation  - 
Part  2;  Laminated  Plates,’  ASME  J.  Appl.  Mech.,  44,  669-676  (1977). 

8.  Hildebrand,  F.  B.,  Reissner,  E.,  and  Thomas,  G.  B.,  "Notes  on  the  Foundations  of  the  Theory  of 
Small  Displacements  of  Orthotropic  Shells,’  NACA,  Technical  Note  No.  1833  (1949). 

9.  Reddy,  J.  N.,  ’A  Simple  Higher-Order  Theory  for  Laminated  Composite  Plates,’  ASME  J.  Appl. 
Mech.,  5,  745-752  (1984). 

10.  Reissner,  E.,  ’On  a  Certain  Mixed  Variational  Theorem  and  a  Proposed  Application,”  Int.  J.  Num. 
Meth.  Engrg.,  20,  1366-1368  (1984). 

11.  Murakami,  H.,  ’A  Laminated  Composite  Plate  Theory  with  Improved  In-Plane  Responses,"  sub¬ 
mitted  for  publication  (1985). 


C4 


- 1  C,‘f ' 


,  1  ,  I>4 -  i  ^11 '  (-  D*  Kt-  +  4 

qU  a 

(A4 

>(*)  -(*)s  , 

I"»  ”  4.  ^  »  <*>J  »<*»  4-  «  (*>S  »(*) 

16  ^6"“  "  " 

“  I  Cu  (t)7  ■!  c 

44g  +  16  «o  n  +  4  no  "  +  «o  " 
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•  Matrix  and  Vectors  V^, ...,  in  Eq.  (21): 
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•  Matrices  [A\], ...,  IC3I  in  Eqs.  (26)  and  (27): 
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•  Vectors  Xi . in  Eqs.  (26)  and  (27): 
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The  component  of  the  vectors  Oi,  ...,  S3  appearing  in  (A8)  are  given  by 
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Table  1  Central  Deflection  uj  for  Symmetric  Cross-Ply  Laminates  in  Cylindrical  Bending  Under 
Sinusoidal  Loading 


Table  2  Central  Deflection  uj  for  Asymmetric  Cross-Ply  Laminates  in  Cylindrical 


Bending  Under  Sinusoidal  Loading 
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Exact  Solution  [I] 
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Present  Theory 
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2.181 

First-Order  Zig-2^g  1 10] 
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1.934 

LCW(71 
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1.979 

-  28  - 


FIGURE  LEGENDS 


FIGURE  1;  Plate  geometry,  coordinate  system  and  trial  in¬ 

plane  displacements. 

FIGURE  2a:  Thickness  variation  of  in-plane  displacement 

of  a  symmetric  3-layer  cross-ply  laminate  for 

S  =  4. 

FIGURE  2b:  Thickness  variation  of  normal  stress  of  a 

symmetric  3-layer  cross  ply  laminate  for  S  =  4. 

FIGURE  3a:  Thickness  variation  of  in-plane  displacement 

of  a  symmetric  5-layer  cross-ply  laminate  for"^ 

S  =  4. 

FIGURE  3b:  Thickness  variation  of  normal  stress  oii  of  a 

symmetric  5-layer  cross-ply  laminate  for  S  =  4. 

FIGURE  4a:  Thickness  variation  of  in-plane  displacement 

of  a  symmetric  5-layer  cross-ply  laminate  for 
S  =  6. 

FIGURE  4b:  Thickness  variation  of  normal  stress  of  a 

symmetric  5-layer  cross-ply  laminate  for  s  =  6. 

FIGURE  5a:  Thickness  variation  of  in-plane  displacement 

of  a  symmetric  9-layer  cross-ply  laminate  for"^ 

>S  *=  4. 

FIGURE  5b:  Thickness  variation  of  normal  stress  of  a 

symmetric  9-layer  cross-ply  laminate  for  5=4. 

FIGURE  6a:  Thickness  variation  of  in-plane  displacement 

of  a  symmetric  9-layer  cross-ply  laminate  for'*’ 

S  =  6. 

FIGURE  6b:  Thickness  variation  of  normal  stress  of  a 

symmetric  9-layer  cross-ply  laminate  for  »  =6. 

FIGURE  7a;  Thickness  variation  of  in-plane  displacement 

of  an  asymmetric  4-layer  cross-ply  laminate  for 

S  =  4. 

Thickness  variation  of  normal  stress  of 

asymmetric  4-layer  cross-ply  laminate  for  S  =  4. 


FIGURE  7b; 


FIGURE  8a: 


Thickness  variation  of  in-plane  displacement 
of  an  asymmetric  4-layer  cross-ply  laminate  for 
S  =  6. 

FIGURE  8b:  Thickness  variation  of  normal  stress  of  an 

asymmetric  4-layer  cross-ply  laminate  for  S  =  6. 

FIGURE  9a:  Thickness  variation  of  in-plane  displacement 

of  an  asymmetric  8-layer  cross-ply  laminate  for  S 

=  4. 

Thickness  variation  of  normal  stress  of  an 

asymmetric  8-layer  cross-ply  laminate  for  S  =  4. 


FIGURE  9b: 
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